Full counting statistics for noninteracting fermions: Exact results and the 

Levitov-Lesovik formula 
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Exact numerical results for the full counting statistics (FCS) for a one-dimensional tight-binding 
model of noninteracting electrons are presented without using an idealized measuring device. The 
two initially separate subsystems are connected at t = and the exact time evolution for the large 
but finite combined system is obtained numerically via a finite dimensional determinant. Even 
for surprisingly short times the approximate description of the time evolution with the help of 
scattering states agrees well with the exact result for the local current matrix elements and the 
FCS. An additional approximation has to be made to recover the Levitov-Lesovik formula in the 
limit where the system size becomes infinite and afterwards the long time limit is addressed. The new 
derivation of the Levitov-Lesovik formula is generalized to more general geometries like a Y-junction 
enclosing a magnetic flux. 



I. INTRODUCTION 

The theory of noise in quantum transport in meso- 
scopic systems is a very active field of research^. In ad- 
dition to the first few moments of the transmitted charge 
the full probability distribution can be studied, called 
full counting statistics (FCS). Neglecting the electron- 
electron interaction Levitov and Lesovik presented an an- 
alytical result in the long time limits using the cumulant 
generating function. It was later questioned pointing out 
that different measurement procedures may lead to dif- 
ferent answers for the distribution function 5 . Recently 
also results for models including the e-e-interaction lo- 
cally were obtained^. 

For the case of perfect transmission the Levitov- 
Lesovik formula-^ yields at zero temperature a delta 
function for the probability distribution w(t, Q) of the 
transmitted charge with the position increasing linearly 
with time. This is in obvious contradiction to the fact 
that w should only have weight at integer values. In or- 
der to obtain the correct answer also in this limit it is 
desirable to start from a formally exact expression. Such 
a starting point was presented by Klicb^ with the help 
of a formula to calculate the trace over Fock space for 
products of exponentials of one-particle operators. 

Klich's formula can be viewed as the general expression 
for the counting statistics of noninteracting fermions, at 
any given time, and without any approximation^. Nev- 
ertheless Klich's attempt to derive the Levitov-Lesovik 
formula from his Eq. (13) is not really satisfying. He 
performes the infinite time limit to introduce the scat- 
tering matrix and notices that therefore the information 
about the length of the time interval has to be put in 
by hand^. As an alternative Muzykantskii and Adamov 9 
start from Klich's formula and use results from the the- 
ory of singular integral equations to proceed. They re- 
strict themselves to zero temperature and switch on the 
bias adiabatically. For the special case of an energy inde- 
pendent transmission probability they recover the T = 
version of the Levitov-Lesovik formula from the solution 



of a matrix Riemann-Hilbert problem. From the approx- 
imate solution of a more difficult Riemann-Hilbert prob- 
lem they also obtain the subleading correction for their 
special case. An extension of this approach to finite tem- 
peratures was presented by Braunecker^ 

The present paper takes a new look at the problem in 
several ways. We first describe the "experimental" setup 
and its formal theoretical description at zero and finite 
temperatures. For noninteracting fermions the central 
object is the time dependent projection operator Pn(t) 
onto one of the subsystems in the Hilbert space of one 
particle. By numerically calculating Pn(t) in section 
III we present exact results for w(t, Q) at T = for a 
one-dimensional tight-binding model. For perfect trans- 
mission the distribution has a nonzero width increasing 
logarithmically with time. In section IV we introduce 
scattering states to describe the long time limit. By us- 
ing the current operator at the connection of the sub- 
sytems we obtain an excellent approximation for Pn(t). 
The new derivation of the Levitov-Lesovik formula from 
Klich's exact formal expression for arbitrary tempera- 
tures presented in section V clearly shows the approx- 
imations involved. We also generalize the derivation to 
more complex geometries like a Y-junction which encloses 
a magnetic flux. In order to make the paper selfcontaincd 
we present an alternative proof of Klich's formula using 
Wick's theorem in Appendix A. 



II. COUNTING STATISTICS FOR 
NONINTERACTING ELECTRONS 

A. Zero temperature limit 

In the following we consider a system which consists of 
two initially separate subsystems described by the Hamil- 
tonians i?o,a with a = L, R. The labels stand for the 
"left" and "right" one-dimensional subsystems treated 
later. The number of electrons in the subsystems in the 
initial state are Nq^ and the total number is given by 
Ntot = Nq,l + No,r. We assume the intial state |$(0)) 
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to be an eigenstate of the i?o,a 

\m) = \e»^)®\e»°<k) . (i) 

The time evolution for times greater than zero is de- 
scribed by the Hamiltonian 

H = H 0<L + H 0tR + V LR = H + Vlr ■ (2) 

The term Vlr which couples the two subsystems will 
be specified later. The probability distribution that Q 
electrons are transfered to the right system after time t 
is given by 

w R (t,Q) = mt)\S[Q-(Af R -N 0iR )]mt)) 

= (*(0)|*[Q - (A/*(i) - JVo,.r)]|$(0)>. (3) 

Here Af R ist the particle number operator of the right 
system and Af R (t) is the corresponding operator in the 
Heisenberg picture. This probability distribution can be 
measured (conceptually simple by "weighing" ) when the 
connection between the subsystems is cut again at time 
t. 

The integral represention of the delta function yields a 
form more convenient for the actual calculation 

w R (t, Q) = ^J d\e- a Qg R (t, A) , (4) 

where the information about the FCS is encoded in the 
characteristic function 

g R (t,X) = ($(0)|e lA ^ ( * ) |$(0))e- 4AWo ' fi . (5) 

This result is valid also for interacting electrons. 

The evaluation simplifies considerably for noninteract- 
ing electrons as the expectation value in Eq. ([5]) can be 
expressed as a iVtot x N to t determinant in the space of the 
initially occupied one-particle states.— As Af R (t) is a one- 
particle operator e lAA ^(*) acts on the one-particle states 
occupied in the initial Slater determinant as e lXPa ^ 
where P R {t) is the time dependent projection operator 

P R (t) = e iht P R e- m = P 2 R {t) (6) 

with P R the projection operator in the Hilbert space of 
one particle onto the states in the right reservoir 

Pr = E \ k f > =Z>> i *> K R \ (7) 

j m 

and h ist the full Hamiltonian in the Hilbert space of 
a single particle. The \k p ) describe standing waves in 
the right subsystem and the m in \m, R) denote the site 
indices. 

If we introduce the projection operator no 
no = X] \ ka )fk,a(k a \ = n ,L + n 0tR , 

k,a 



where the fk, a are the (T — 0) Fermi functions of the 
subsystems in the initial state the operator 

P R (t) = n P R (t)n Q . (8) 

determines the zero temperature FCS. Using e lAp n(*) 
1 + ( e iA - l)P R (t) we obtain 

g R (t,X) = e-^^detn, [n + (e tX - l)P R (t)] ,(9) 

where detn is a determinant in the subspace of the ini- 
tially occupied one particle states. If we denote the eigen- 
values of P R (t) by p m (t) the characteristic function takes 
the form 

JV tot 

g R (t,X) = e" aJV °.* J] [l + (e^-l) Pm (t)} 

m— 1 

JVtot 

= ^ c n (t)e^ n - No ^ . (10) 

n=0 

The coefficients c n (t) can be calculated recursively from 
the product representation. This leads to the probability 
distribution 

Ntot 

w R (t, Q) = Y, c ^) 5 iQ - ( n - n o,r)) ■ ( n ) 

To obtain exact results for the FCS one first has to cal- 
culate P R (t) and then obtain its eigenvalues p m (t) . 

B. Finite temperatures 

The starting point for finite temperatures is again 
given by Eq. ([5|). We take the factor e~ tXN °- R inside 
the expectation value and replace it by e - 4A/ H(*)\ \y c 
then perform a grand canonical averaging corresponding 
to the statistical operator 

g-/3i(Ho,i,-/i£A/i,) g— Pr(H 0iR — ixrNr) 

= e-"°/Z (12) 
The probability distribution w is then given by 

t«£(f,Q) = i- / d\e- lX Q (e'^C'e-'*) , (13) 

where (•) denotes the averaging with the statistical op- 
erator pq. The information about the FCS is therefore 
encoded in the characteristic function 

g R (t,X) = (e lXArR{t) e'^) . (14) 

Again this result is valid also for interacting electrons. 

For noninteracting fermions the characteristic function 
can again be expressed as a determinant but now in the 
full one particle Hilbert space. One can either use Klich's 
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trace formula^ (Eq. (|A13[) ) or use the proof using Wick's 
theorem presented in apppendix A. This yields 



gn(t, A) = det 



1 



1 n Q 



(15) 



where the Fermi functions in the definition of no are now 
the ones for finite temperatures. This implies that h 
is no longer a projection operator. At T — Eq. (fK 
reduces to Eq. ©. 

As the next step we rewrite Eq. (|15|) in the form 



Q R (t, A) = det 



1 



n 



, e iXP R {t) /= + e -*A= ) 



In contrast to statements in the literature^ it is useful to 
write P R (t) as a sum 

Pn{t) = Pr(0) + I P R {t')dt' = P R + 6P R (t) . (16) 



With e lAP «(*) = 1 + ( e ' A - l)(P R + SP R (t)) we obtain 



lXPR(t) ( n , L + e~ lX h QM ) - no 

= 5P R (t) [no,L {e lX - 1) - 
= SP R (t)[a L (X)-a R (X)] 



(17) 



«o,i?. (e 



1)] 



The general result for the characteristic function then 
takes the form 

9 R (t, A) = det [l + SP R (t) (a L (\) - a R (\))] . (18) 

As shown in section VI this is the starting point for 
a transparent derivation of the approximate (long time 
limit) Levitov-Lesovik formula 

= ±J de\n{l + T( £ )[(e iA - 1)/ L (e)(l - f R (e)) 

+ (e- tX -l)f R (e)(l-f L (e))}}, (19) 

where T(e) is the transmission probability for the scat- 
tering of a single particle between the two subsystems. 

Note that compared to the zero temperature result the 
operator added to the unit operator in the determinant 
in Eq. (j 1 8|) cannot be written as a product of a func- 
tion of A and an operator independent of A. Therefore 
the numerical effort for an exact numerical calculation of 
w R (t, Q) is much higher than at T = 0. 



The unperturbed one particle Hamiltonians of the sub- 
systems are given by 



2J T m (\m)(m+ 1| +H.c.) + v 



m=-(Mi-l) 
n Mr-1 



ho, R = -j J2 (\m){m + l\ + H.c.) . 



(20) 



The number of sites m in the subsystems are given 
by M a . In the explicit calculations we later special- 
ize to r m = B/4 for m = — (Ml — 1), —2 and r_i = t^. 

The one particle eigenstates |A^) of the right subsys- 
tem are standing waves 



(m\kf) = 



1 



J* 



1 



(21) 



with the integer j e [1, M R \. The corresponding energies 
are given by e k R = — (B/2) cos fc p , i.e. the total band 
width is given by B. 

In order to obtain a stationary current flow in the long 
time limit t — > 00 one first has to take the limit M a — * 00. 
In this limit the quantum numbers kj become continuous 
with k E [0, 7r] and one has to go over to the delta function 
normalization of the states 



Mr + 1 



k; 



\k,R) 



(22) 



with (k,R\k',R) =5{k-k'). 

In the left subsystem the eigenstates \kf) are pure sine- 
waves only in the special case i>o = and r m = const. 

The coupling between the two subsystems is described 
by the hopping term 



v LR = -t R (\0)(l\+H.c) 



(23) 



The special case i>o = and r m = const, corresponds to 
two subsystems which, apart from their initial filling, are 
identical. The numerical results presented in this section 
are for this simplest case. The model with r_i = ti 
and taking arbitrary values is a model which allows 
resonant scattering through the "quantum dot" at site 
zero. Some properties of this model simplify in the "wide 
band limit" B — > 00. 



III. EXACT NUMERICAL RESULTS FOR T — 

A. The model 

In this section we present exact numerical results for 
the probability distribution w(t, Q) for a one dimen- 
sional tight binding model. The generalization for the 
model with more than two leads discussed in section VI 
is straightforward. 



B. Results for finite times 

As discussed in the previous section one has to address 
the eigenvalue problem of the operator P R (t) in order to 
calculate the probability distribution w(t,Q). For the 
calculation of the eigenvalues p m (t) of P R (t) we use 

M R 

(WP R (t)\4) = 5>£(t)MH^W> , (24) 

771=1 



20 40 60 

m 

FIG. 1: Eigenvalues p m (t) of Pn{t) for a system with Ml = 
Mr = 100, N ,l = 80,iVo,H = 50 and t = tB/A = 40 for 
tR = AtR/B = 1.0 (circles) and tR — 0.5 (squares). Starting 
with m = 51 they are compared to the NdiB eigenvalues of 
PR,dift(t) defined in the text: tR = 1.0 (diamonds) and tR — 
0.5 (triangles up). 



Q 

FIG. 2: Weights of the probability ditribution WR(t,Q) for 
the two parameter sets of Fig. [1] From the eigenvalues of 
Pit{t): tR = 0.5 (filled squares), tR — 1.0 (filled circles) and 
the eigenvalues of pR >( us(t): tR = 0.5 (triangles up), tR = 1.0 
(diamonds). 



where p,, p! label the occupied one particle states. With 
the help of the spectral representation of the full Hamil- 
tonian h 

h = ^2\ki)e ki {ki\ (25) 

the time dependent quantities in Eq. (|2~4"|) can be ex- 
pressed as sums 

{m\k a ;(t))=Y^{rn\k i )e- i ^ t {k i \k a ;) . (26) 

i 

In order to obtain the eigenvalues the resulting 
iVtot x ^Vtot matrix has to be diagonalized numerically. 

As discussed in more detail in the next section the 
results are almost independent of the size of the system 
and only depend on the values of the chemical potential 
Hl and piR for times smaller than the time it takes the 
"charge fronts" which move into the subsystems after 
connecting them to return to the the connection point af- 
ter having been reflected at the boundaries of the system. 

We always consider the case ATq.l > No t R, i.e. the net 
particle flow is from left to right. As Pn(t) is a projection 
operator the eigenvalues p m {t) of Pn(t) obey the bounds 
< p m (t) < 1. In the figures we show the eigenvalues 
p m {t) in descending order. 

For t = the matrix Pr(0) is diagonal with Nq^r eigen- 
values 1 and Nq >l eigenvalues 0. For t > the Nq^r 
eigenvalues almost identical to 1 persist. As shown in 
Fig. [1] for perfect transmission N t additional eigenvalues 
almost identical to 1 occur where the value N t ~ t is dis- 
cussed in detail laterJ^ The crossover from eigenvalues 
close to one to those close to zero is smooth. For non- 
perfect transmission the eigenvalues which significantly 



differ from 1 and start at Nr^ + 1 and approximately 
extend to Nr^ + N t . 

If one has a look at the Levitov-Lesovik formula Eq. 
(0 for T = 

\ng { B L \t,\) = l- deln[l + T(e)(e^-l)] . (27) 

one expects the states below [ir which are initially occu- 
pied in both subsystems to be irrelevant. We therefore 
also calculated the eigenvalues of the N^m x iVdig sub- 
matrix PR,difi(t) in the space of the Ndis = Nq^ l — No,r 
states of the left system in the energy window [/ir, p L }. 
As can be seen in Fig. [1] these eigenvalues deviate little 
from the eigenvalues of the full N to t x AT tot matrix. 

The weights of the probability distribution correspond- 
ing to the parameters of Fig. Q]are shown in Fig. [2] . For 
the non-ideal transmission (tR = Hr/B = 0.5) the distri- 
bution calculated from the "extra electrons" only, agrees 
very well with the one from the full calculation. For the 
perfect transmission case slight deviations show up. 

If one doubles the size of the system and the number of 
electrons but leaves t fixed, the result for WR(t, Q) is al- 
most identical to the one shown in Fig. [2] Our approach 
using finite systems and nevertheless see "asymtotic" 
properties in a time window is very useful also for in- 
teracting fermions for which effective methods have been 
worked out recently to describe the time dependence for 
finite systems. Promising techniques are the time depen- 
dent density matrix renormalization group (DMRG) 13 ' 14 
and the time dependent numerical renormalization group 
(NRG)^ 



Note that the T = Levitov-Lesovik formula predicts 
zero width of the probability distribution for the perfect 
transmission case which is in contrast to the exact nu- 
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FIG. 3: Second order cumulant K2 for perfect transmission at 
T = as a function of time for the system parameters of Figs. 
Q] and [2] The circles show the result using the matrix Pa(t), 
while the squares are the results for the matrix PR,dis(f) which 
only treats the electrons in the energy window [fiR, 



merical result in Fig. [21 To study the non-zero width 
as a function of time one can calculate the second order 
cumulant. 



C. Low order cumulants 

The low order cumulants Kj of wn(t,Q) are obtained 
from the Taylor expansion of In gn(t, Q). Using det a = 
exp (tr In a) in Eq. ^ one obtains for i up to 4 





= trP R 


K 2 


= ixPn 


K3 


= ixP R 


«4 


= trP R 



K 

trPi 
StrPl 



(28) 



2trP| 
12trP| 



6trP, 



R ' 



In the long time limit discussed in the next section an- 
alytical results for the low order cumulants can be ob- 
tained. 

In Fig. [3]we show results for k 2 for perfect transmission 
and the parameters of Fig. Q] . The circles correspond 
to the full calculation and the squares to the result using 
-Pfl.diff (t). It will be shown analytically in the next section 
that K2 increases logarithmically at T = in the long 
time limit for perfect transmission. The numerical results 
for Ki with i > 3 show no logarithmic increase. Apart 
from small transient effects these cumulants stay close to 
zero. 

The Levitov-Lesovik formula for T = predicts zero 
width of the probability distribution also for the zero 
bias case (II — M-R f° r arbitrary transmission. Absolutely 
no charge transport is predicted in contrast to the exact 
numerical results not shown here. 



IV. THE LONG TIME LIMIT 
A. Introduction of scattering states 

We now examine if Pr(<) can be discussed analytically 
in the long time limit. For the calculation of the matrix 
elements of Pn(t) it is useful to split off a phase factor in 
the time dependent one particle states and dehne 



|**> t = e fc V|tf(i)> 



-'(^-eL)* 



(29) 



It is known from scattering theory- that in the long time 
limit (after taking the limit M a — > oo) the states \k,a)t 
converge to the scattering states. 
This can be seen using an Abelian limit procedure 



lime %{h th - a)t \k,a) = lim 77 / < 

t~>oo r f ^0 J Q 

VI] 



■•»*e- i(fc - 6 *-- )t |*r ) a)<ft 



lim 

\k,a+) 



(30) 



The state on the rhs of the last equality is just the scat- 
tering state with outgoing scattered waves. 

The explicit form of the scattering states for our Hamil- 
tonian h with r m = B/4 for m < —1 is given by 



(m\k, L- 



1 



/2tt 
f 



'2tt 



j. Ahra 

t k e 



m < — 1 



m > 1 



(31) 



with rfc and tk the reflection and transmission amplitude. 
Similarly the scattering states with incoming wave from 
the right are given by 



(m\k,R+) = 



1 



/27T 
1 



-ikm _|_ f ke ikm 



/2tt 



-tke 



-ikm 



m > 1 
m < -1 . (32) 



If one returns to the description of the hnite system at 
finite times the question arises how well the approxima- 
tion 



\k](t)} 



N a + 1 



\kj, a+)( 



(33) 



works. 



In Fig. [4] we show 2(N L + l)\(m\kf) t \ 2 for different 
values to > of the site index for identical subsystems 
with Ml = Mr = 100, hopping matrix element tR = 
^r/B = 0.5 and a state in the middle of the band. In the 
approximation of Eq. (|33p this quantity should be time 
independent and given by the transmission probability 
ItfcJ 2 . The figure shows that the larger the value of m 
the later the oscillatory convergence sets in. 



6 




15 20 
tfl/4 



25 30 



FIG. 4: This figure shows how the long time approximation 
in Eq. ([33|) works for finite times. The quantity plotted should 
converge to the value \t kj | 2 for all values m > of the site 
index. Shown is the result for a state in the middle of the 
band for a system with Ml = Mr — 100 and tn = Atn/ B = 
0.5. For the times shown the results agree within the plotting 
accuracy with corresponding calculations for Ml = Mr 
100. 



useful approximation also for finite systems 



(k°\P R (t)\kf,) « S^S^rSjj, 

(kj,a+ \ja^i\kj',a'+) > 



(37) 



^(M + l){M a , + 1) i(e*j,.-«v,a') 

For the calculation of the current matrix elements we 
also have to know how (0|fc, a+) is related to the scatter- 
ing amplitudes in Eqs. (|67|) and (|32|) . The Schrodinger 
equation (l|/i|fc,a+) = efc )0 (l|fc, a+) yields 



t R (Q\k,L- 



tk B t R (0\k,R+)= 1 -±^? r . (38) 



/2tt 4 



/27T 4 



This allows to express the current matrix elements in 
Eq.([3"7p in terms of the t k and r k , e.g. 

(fe,i+|io-i|fc',i+) = ^t* k t k ,(e-* k -e* k ')^ . (39) 



and 



It would be totally wrong to use Eq. (|33[) in Eq. 
for all sites of the subsystems. The only time dependence 
would then result from the factors e lek »> at e te V,a' . 
They could be gauged away by absorbing them in the 
definition of the basis states making the eigenvalues 
p m (t) time independent. This would imply that g R (t, A) 
and with it the FCS is time independent. This shows 
clearly that one should apply Eq. (|33|) only to the 
calculation of observables localized in the neighborhood 
of the junction. 

It is therefore useful to express P R (t) as an integral 
over the current operator between the two subsystems as 
in Eq. (H 



P R (t) = Pr(0) + [ j^i{t')dt' = P R + 6P R (t) (34) 
Jo 

with the current operator given by 

jo^i = it* (|1> <0| - |0)<1|) . (35) 
The matrix elements of P R (t) are therefore given by 



(k,L+\ J0 ^ 1 \k\R+) = —[ t * k r k ,(e- ik -e ik ') 

Zir 



+ t* k (e 



- ik 



e~ lk '))j- (40) 



These expressions simplify in the in the wide band limit 
— fJ, R \ <C B for the scattering states in the energy 
window between /zl and /j, r . Then one can replace k 
and k' by JtF = (&f,l + fcF,ii)/2. This yields 



(k,L+\j a ^i\k',L+) « V -^t* k t k , 



(k,L + \j ^i\k',R- 



2_7T 

2tt 



J.* ~ 



(41) 



with the Fermi velocity vf = v(kp) = 2BsmkF- As 
the transmission amplitudes may have a rapid variation 
in this energy window the replacement k, k' — ► kF has 
not been made for the factors involving the scattering 
amplitudes. 



B. Low order cumulants 



(k]\P R (t)\k],) = S a ,a'Sa,RSjj' 

W{t')\3 Q ^i\4(t'))dt' . (36) 



As expected from Fig. Q]the use of the approximation 
in Eq. (|33[) to evaluate the matrix elements of the local 
observable Jo— >i works surprisingly well even for rather 
short times. With this replacement the time integration 
in Eq. (|36|) can be carried out and one obtains the very 



As for the Hamiltonian in Eq. ([20|) the transmission and 
reflection amplitudes can easily be calculated (see the ap- 
pendix), the matrix elements of P R (t) in the scattering 
approximation Eq. (|37|) . called P R (t) in the following, 
are known explicitely. Without further approximations 
the eigenvalues of this matrix have to be calculated nu- 
merically. The low order cumulants /tj of w R (t, Q) on the 
other hand can be discussed analytically. We begin with 
the T — limit where the cumulants have the simple 
form presented in Eq. 
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Results for T — 



Using Eqs. ([28)1 and (j37|) the explicit form of the time 
derivative of the first order cumulant for T = is given 
by 

~ H]^7-Y( fc M^+ljo^i|^,a+> (42) 

For symmetric subsystems (r m = -B/4, v = and 
Mx, = Mr) the contributions of the states with energy 
smaller than hr{< (j-l) which are occupied in both sub- 
systems exactly cancel and only the states \kj^) in the 
energy window \pr, fii] contribute. In the thermody- 
namic limit n/{N a + l)/£ fco 0) -> f(-)dk = J(-)de/v{e) 
the cancellation of the contributions from the states be- 
low hr holds generally in the approximation of using the 
scattering states and one obtains using Eq 



1 f^ L 
Al t "5-/ T(e)de. 



(43) 



This is the zero temperature version of the Landauer- 



= t 



fc(e)| 



the transmission 



Biittiker formula with T(e) 
probability*^ 

In Fig. [5] we show how the long time result for ki 
is approached at finite times. From the results in Fig. 
[4] one expects also an oscillatory approach of the exact 
numerical result to the Landauer-Buttiker formula. This 
is confirmed in Fig. for different values of the hopping 
parameter between the two subsystems. 
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FIG. 5: Comparison of the exact and the approximate long 
time result for d,K\/dt given in Eq. (|42[1 for a system with 
Ml = Mr = 200, N ,l = 120, Nq,r = 100 and three different 
values of ta. For the smaller values of tu the convergence to 
the approximate constant result (dotted curves) is faster. 

For perfect transmission the long time approximation 
for the first cumulant is given by 

ki - (44) 



with A/x = fiL — M-R- Neglecting the transition region in 
Fig. 1 this is realized by N t eigenvalues 1 of -P^diff which 
implies 



Ni 



t 

2^ 



(45) 



The numerical evidence that the number of eigenvalues 
Pm,diff(^) which are different from zero is in general ap- 
proximately given by Nt yields for integer N t 

N diB N t 

K\ = Pm,diff(t) ~ 2J Pm,6AB{t) ■ (46) 

m— 1 m— 1 

This can be compared to the approximate evaluation of 
Eq. (|43| using the trapezoidal rule for Nt subintervals. 
Together with Eq. (|45| one obtains 

N t -1/2 Nt 

m— 1 ' m— 1 

It is therefore tempting to conclude 

m — 1/2 

Pm.diff(t) « T(fi R + t^A/z) (48) 

for m — 1, Nt and zero otherwise. With Eq. (|4H>1) the 
argument of T in this equation can also be written as 
[Xr + 2n(rn — 1/2) /t. To really justify the guess in Eq. 
(|48p one has to examine also the higher order cumulants 
or compare with the exact numerical results. 

For T — the explicit form of the second order cumu- 
lant is 

K 2 = T,T,K\ p n(t)\K)(K\PR(t)K) 

/i, a a, a' 

.2 

-\(k a ,a' + |j ^i|^> a +)l 

kfj, ,a k a ,a 

'sin [(e k 



^ ^, (JV + l)(jV , + l)' 



(c fc -e fc J/2 J 



(49) 



where fi labels occupied and a unoccupied one particle 
states. 

If one performs the thermodynamic limit and replaces 
the last factor on the rhs of Eq. (|4l?l) by 2irt5(ek fi — efc Q ) 
and uses Eq. (|40|) one obtains Lesovik's shot noise result^ 



^2 



t 

2^ 



T(e)(l-T(e))de 



(50) 



which is the result for ki which holds for the zero tem- 
perature Levitov-Lesovik distribution in Eq. (|27|). It is 
consistent with the guess in Eq. (|4"8")) . To further test 
this approximation for the eigenvalues p m ,dis{t) we show 
a comparison with exact results in Fig. [6] Note that 
the aproximation p m ,dis(t) ~ T m {t) fails to to describe 
the transition region to the zero eigenvalues. This is in 
contrast to the eigenvalues Pmdiff(t) usm g the scattering 
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m 



FIG. 6: Eigenvalues PR,dm(t) for a system with Ml = Mr = 
200, N , L = 140, N ,r = 100 and t R = 0.3 at t = 80 (filled 
circles ), t — 120 (filled diamonds ), and t = 160 (filled trian- 
gles ). For t = 80 the eigenvalues pit difi(i) ar e shown as open 
squares and for t = 160 the T m defined in Eq. (|47|) as open 
circles. 



approximation Eq. ([57)1 which at t = 80 show an overall 
good agreement with the exact p m ,dis{t)- 

For the special case of perfect transmission T(e) = 1 
the approximate result in Eq. (|50j) yields zero width for 
the probability distribution wn(t,Q) in contrast to the 
exact numerical result shown in Fig. [TJ 

In order to obtain a non zero result for k 2 for the per- 
fect transmission case at T — the integrations have to 
be performed more carefully. In order to keep the cal- 
culation as simple as possible we use the approximation 
to only treat the electrons in the energy window [/j,r, fit] 
which coresponds to Nq,r — 0. In addition we work in 
the wide band limit and assume an energy independent 
transmission probability T, i.e. approximate the current 
matrix element (k,L+ \jo-ti\k' , L+) by Tvp/(2ir) (see 
Eq. (HI])). This yields for the trace of P R (t) entering the 
expression for K2 in Eq. I]28p 

2T 2 f tAf * /2 sin 2 u , 
-> TtrP H (t) 5- / . (51) 

The remaining integral is logarithmically divergent in the 
long time limit t ^> 1/B and one obtains 

tAu - - f 2 
^2 w — ^T(l - T) + — ln(tA///2) . (52) 

This agrees with the corresponding result from the ap- 
proximate solution of the matrix Riemann-Hilbert prob- 
lem in Ref. [H. The logarithmic correction term is only 
important for the case of (nearly) perfect transmission. 



As already mentioned for the numerical result using the 
exact p m {t) the cumulants ftj for i > 3 vanish (apart 
from small transients) for perfect transmission. There- 
fore the values of WR(Q,t) lie on a Gaussian of width 
~ [ln(tA/i/2)] 1 / 2 to an excellent approximation. This is 
in agreement with the general result for the sublcading 
correction for the special case considered in Ref. |{§. 

2. Finite temperature results 

As the finite temperature result for gR(t, A) in Eq. (fl~8)) 
contains two different functions of A the calculation of the 
cumulants is more tedious than for T = 0. 

/si(i) = tt[6P R (t)no] (53) 
K2 (t) = tx[5P R (t)5n- 5P R (t)noSP R (t)no} 

= tr [5 PR(t)(i-n )5P R (t)n Q ] (54) 

with Sn = fiQ >L — no t R. In the second equality for «2 
we have used the fact that P R (t) = Pr + SP R (t) is a 
projection operator which reads for 5PR(t) 

SP R (t) = P R SP R (t) + SP R (t)P R + 8P 2 R {t) . (55) 

The long time result for k,\ is just the Landauer-Buttiker 
formula integrated over time 

«! = £<fc?|M>fl(t)|fc?>/£ (56) 

^ hi T (t)(fL(e)~fR(e))de. 
The second form for K2 shows explicitely its positivity 
^ = T,T,\W\ SP *( t M')\ 2 fk i ( 1 - & (57) 

i,a j,a' 

* hi de{T2{e) [/i(e)(1 " /l(£)) + /fl(e)(1 ~ /fl(e))] 

+ T(e)(l - T(e)) ^(0(1 - + / fl (c)(l - f L (e))}} 

The energy integral is the contribution linear in time and 
is obtained by the generalization of the steps to derive Eq. 
(|50p . At finite temperatures there is a contribution linear 
in t even for the perfect transmission case, which masks 
the logarithmic contribution in Eq. ([52| for long enough 
times. 

We show in the next section that the Levitov-Lesovik 
approximation Eq. (|19[) amounts to include the linear in 
t contributions in all orders. 



V. DERIVATION OF THE LEVITOV-LESOVIK 
FORMULA 

A. The one-dimensional model 

After we have shown how the use of scattering states 
simplifies the theoretical description in the long time 
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limit we can use Eq. (|18[) to discuss which approxima- 
tions are involved to derive the Levitov-Lesovik formula 
Eq. (fT9|) . If we define the operator b in the single particle 
Hilbert space 

b(t,X) = SP R (t)[a L (X)-a R (X)] (58) 

the exact result for the characteristic function reads 



m-\-l 



g R (t, A) = e 




tr b n 



(59) 



In the following we simplify the exact expression for tr b m 
in two steps. 

In the first step we introduce the scattering states using 
Eq. ([57)) . In the second step after taking the thermo- 
dynamic limit the product of the factors (e zej,i * — 1)/ 
with e 3: ; = ej — ei is used to obtain a product of "energy 
conserving" delta functions 



1 e «i 2 ,i 3 f 



f 



ic. 



IC, 



(60) 



sin (e il<i2 t/2) sin (e i2l i 3 t/2) sin (e imiil t/2) 



til ,il 



,/2 



e Wl/ 2 



-> (27r) m ^(e^ - e i2 )(5(e l2 - e l3 ) . . .6(e im - e n ) 

The k- integrations in the thermodynamic limit can be 
converted to energy integrations. Because of the product 
of delta functions only one energy integration remains. 
The summations over the quantum numbers L, R still 
have to be performed. They can be expressed as a trace 
in a 2-dimcnsional space spanned by a = L, R 



tr 6" 



t 

2^ 



de tr (2) [c( 2 )(e, A)] r 



(61) 



with 



c (2) (e,A)= (62) 

2tt f(k,L+ \j\k, L+)a L -(k,L+ \j\k, R+)a R 
v (e) \(k, R + \j\k, L+)a L -(k,R+ \j\k, R+)a R/ 

= / \t k \ 2 a L (e,X) -t* k r k a R (e,\T 
\t k f%a L (e,X) \t k \ 2 a R (e,X) 

Here k = k(e), a L{R) (e,X) = f L(R) (e)(e ±iA - 1) and we 
have dropped the index — > 1 of the current operator. 
From Eqs. ([ST))) and (|6"Tj) one obtains 



tr[ln(l + 6)] - j-Jdt 



\m+l 



(—'\Y 

tr( 2 )[ C ( 2 )( £ ,A)] f 

rn — 1 



m 



— j rfetr (2) [ln(l 2 + c (2) (e,A))] 



t 

2^ 



deln[det (2) (i 2 +c (2) (e,A))]. (63) 



The 2x2 determinant can easily be calculated. Using 
\r k \ 2 = l-T(e) and a L a R = -(f R a L +f L a R ) one obtains 

det (2) (i 2 + c (2) ) = l + T(e)l(eS x -l)f L (e)(l-f R (e)) 
+ (e-^-l)f R (e)(l-f L (e))} . (64) 



If one inserts this into Eq. (|63|) this completes the deriva- 
tion of the Levitov-Lesovik formula. Note that in the 
derivation it was not necessary to invoke the wide band 
limit. 

If one is interested in logarithmic corrections in time 
to hig R (t, A) the approximation in Eq. (|60j) should not 
be made for carrying out the energy integrations (see the 
discussion of k 2 leading to Eq. (|52)l in section IV b). 



B. Multi-lead geometries 

The FCS in multiterminal circuits was addressed by 
various authors i 19 ' 20 For noninteracting electrons the 
Levitov-Lesovik formula for the quantum wire discussed 
in the previous subsection can easily be generalized to 
such more complex geometries. As an example we show 
in Fig. [Jj a Y-junction pierced by a magnetic flux. The 
zero temperature limit of a Y-junction was already dis- 
cussed in ref. 3. The proof given below works for the 
general case of M legs at finite temperatures. 




1,3 2,3 

-e e- 



FIG. 7: Junction of M = 3 quantum wires connected to a 
"ring" indicated by the filled circles, which is pierced by a 
magnetic flux. In the initial state the hopping matrix elements 
between the ring sites and the adjacent wire sites vanish. 



We denote the states on the central "ring" as |0, a) 
with a — 1, ...M and the states in the legs in the site rep- 
resentation as \m, a) with m = 1, 2, ...oo. The connection 
of the legs with the ring generalizes Eq. (f23|) to 



A I 



J2ra(\0,a)(l,a\+H.c), 



(65) 



where the r a can be assumed to be real. In the legs we 
assume nearest neighbor hopping r = 1. At time t = 
the subsystems with different chemical potentials (and 
temperatures) are connected via the hopping term v in 
Eq. (|65p . Instead of the particle number in the right part 
of the system we now monitor the particle number in leg 
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M. The relevant current operator is then given by 

3o-+i,M = iTM(\l,M)(0,M\ - \0,M)(1,M\) . (66) 
The scattering states take the form 



i,a'\k,a+) = 



1 



/2tt 



-ikm i „ ikm 
+ r a .ke 



) , a = 



/27T 



t 



a' * — a ikm 



, a 7^ a 



(67) 



where m > 1. When the ring encloses a magnetic flux 
\ik generically differs from \ f^ a \ Z^. 

In order to calculate the current matrix elements we 
have to generalize Eq. (|38|) to 



T M (0,M\k,a- 
T M {0,M\k,M- 



1 



2tt 
1 



M<-a 



a^M (68) 



(1 +r M ,fe) 



Having defined the multi-lead model we can now gen- 
eralize our description of the FCS in section III . The 
operator no now takes the form 

M 

no = "o,ring + E no,a = n' + no,M (69) 

and the operator b in the characteristic function 
gM(t, A) = det(l + b(t, A)) is given by 

b(t, A) = 6P M (t) [n' (e iX - 1) - n , M (e~ iX - 1)] . (70) 

As in the preceeding subsection we evaluate the deter- 
minant in the long time limit by calculating tr b ,n to all 
orders. As we are only interested in the linear in t terms 
we can neglect the ring state contribution to the trace. 
Wc then obtain as the generalization of Eq. (|6"Tj) 



tr b m - — 
2tt 



de tr (M) [c (M )(e, A)] r 



(71) 



The current matrix elements entering the matrix 
c (M) (e, A) are given by 



(k,a+ \ jo->i,M\k,a'+) 
(k,a+\j ^ hM \k,M+) 
(k,M+\j ^ M \k ) M+) 



v (k) , ,M<~a\* + M«-o' 

V -P-{t^' a Tr M , k (72) 

Z7T 

«(fc) 



2tt 



(l-\r M , k \ 2 ) 



Here a, a' 7^ M. For the Y-junction with M — 3 the 
matrix I3 + C( 3 ) (e, A) is given by 



l+i*ai tja 2 — t\a^ 

C(3) = I *2 a l 1 + *2 a 2 - *2«3 

r*ai r*a 2 1 + D 



(73) 



Here we have introduced the abbreviations t¥ < ~ a — > t a 
and rjw.fe -> r and defined ai (2 ) = ii(2)/i(2) (e' A - 1), »3 = 



rh{e~ lX - 1) and D = (1 - |r| 2 )/ 3 (e- iA - 1). As for 
general values of M the determinant of this matrix is 
best calculated by expanding in the last row. This yields 

det (13+0(3)) = {l+D){\+t* l ai+t* 2 a2)+r*a s ,{t* l ai+t2a 2 )- 

(74) 

For arbitrary values of M the sum over the t*ai extends 
to M— 1. Therefore the generalization to arbitrary values 
of M reads with d± = e ±lX — 1 and using d+d- = — (d+ + 

Af-l 

det(l M + c (M) ) = 1 + d+(l - / M ) X! I*f "°| 2 /a (75) 

a=l 

+ d_/ M (l - | rjltf)fc | 2 - £ l^f/aj- 



a=l 



The conservation of probability implies 

A/-1 



i-k M , fe i 2 = E i*r 



•A/ 1 2 



(76) 



As generally |<^ f ^ a | 2 7^ |^^ M | 2 it requires an additional 



argument to bring Eq. (|75[) into its final form. It is shown 
in the appendix B that 



A/-1 



Af-l 



(77) 



also holds in the presence of a magnetic flux. Using Eq. 
(|63| for arbitrary values of M we obtain the generalized 



Levitov-Lesovik formula . With T a (e) = \t% r ^ a \ 2 it reads 



with 



In g M (t,X)^ — I deln[l + A(e,A)] (78) 



Af-l 

A(e,X) = E T «( £ ) [(e iA -l)/a(e)(l-/M(e)) 

a=l 

+ (e- a -l)/ M (e)(l-/„(e))] (79) 
Again it was not necessary to invoke the wide band limit. 

VI. SUMMARY 

The aim of this paper was to clarify the validity of 
the Levitov-Lesovik formula for the FCS of noninteract- 
ing fermions and to present a simple derivation which 
clearly shows the aproximations involved. As a starting 
point we used the formula for the characteristic func- 
tion in Eq. (p~5|) first derived by Klich. 8 No measuring 
device is used in this description. As Klich noted his for- 
mal result should be viewed as the general expression for 
the counting statistics of noninteracting fermions, at any 
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given time and without any approximation. Our alterna- 
tive derivation of Klich's formula for finite temperatures 
presented in Appendix A uses Wick's theorem. 

The first step in the actual calculation of the charac- 
teristic function is to obtain an explicit result for the time 
dependent projection operator Pr(£). This was done ex- 
actly for a microscopic lattice model using a numerical 
approach and analytically in the long time limit using 
scattering states. It was shown by a comparison of the 
two aproaches that the analytical expression in Eq. (|5T|) 
works very well after transients with a time scale ~ 1/B 
have died out. At zero temperature the eigenvalues p m {t) 
of the operator Pn(t) directly determine the probability 
distribution w(Q,t). The fact that the Levitov-Lesovik 
formula provides only an approximate description of the 
FCS was studied in detail using low order cumulants. A 
comparison to results obtained by different methods^ was 
given. 

In section VI we presented a straightforward deriva- 
tion of the Levitov-Lesovik formula summing all linear 
in t contributions to \ng R (t, A) using Eqs. (fTo")) and ([57)) . 
An important step in the derivation was to first bring 
Eq. (fTS]) into the form Eq. (|18p introducing the opera- 
tor 6Pfi(t). In the derivation no assumptions about the 
band width an the energy dependence of the transmis- 
sion probability were necessary. The derivation applied 
to the simple "experimental" setup studied in this paper 
shows that the Levitov-Lesovik formula provides an ex- 
cellent approximation in the long time limit, except in the 
case of (near) perfect transmission at very low tempera- 
tures. An extension of this new derivation to multi-lead 
circuits with possible broken time reversal symmetry was 
also given. 
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APPENDIX A: KLICH'S FORMULA USING 
WICK'S THEOREM 

In this appendix we discuss the expection value of an 
exponential of a selfadjoint one-particle operator A 



A 



where the c'P (as well as the c^' ) obey the usual an- 

ticommutation relations {ci,c^} = <Sy,{cf, Cj} = 0. In 
order to avoid mathematical subtleties we assume the di- 
mension of the one-particle Hilbcrt space to be finite (Af). 
The creation operators for the two different one-particle 
basis states obey the usual relation 



W a \3 



<ctc 



E 



a\c\c\ 



(Al) 



1,3 



it) 



(A2) 



We consider expectation values with a statistical oper- 
ator in Fock space 



po = e- Ho /Tre- Ho = e- H °/Z 



with 



Q C„- Ci 



(A3) 



(A4) 



In the section II we use H — Pl{Hq,l — PlNl) + 
/3r(Ho,r — PrNr), where the Hq m are one-particle oper- 
ators and the J\f a are the particle number operators of the 
left and right part of the system. The grand canonical 
part ion function Zq is given by 



Zn = 



det(l 



(A5) 



With the operator in the one particle Hilbert space 

1 ,., 



n 



i 



(A6) 



the expectation values of bilinear operators in the A- 
representation can be written as 



( c W) 



= £<Mli)<c]c J -)(i|A) = ( M |n |A) 



(A7) 



We now consider the expectation value of e aA , where a 
is an arbitrary complex number 

(e aA ) = (e a £* = Y[ ( [l + (e aax - l)c{cx\ ) . 

(A8) 

The product runs over the M different eigenstates labeled 
by A. If one multiplies out the product one has to caclu- 
late expectation values of products of different fi\ = c^ca- 
They are given by Wick's theorem^ as a determinant of 
the bilinear expectation values in Eq. (|A7|) 



(4^1 



=A CO = 



(4i c ai> ■■• (4i c O 
(4 cai> ■■■ (4 m c O 



(A9) 



-1) 



This expectation value is multiplied by Jll^i( e 
These factors can be incorporated into the matrix ele- 
ments on the rhs of Eq. (|A9|) using Eq. (|A7|) 

t 



(e^-l)(cc Aj ) = (A 



l)no|Aj> . (A10) 



A compact expression for (e aA ) is finally obtained by 
comparison with the formula for the determinant of 
1 + 6, where b is an arbitrary linear operator in the M- 
dimcnsional Hilbert space 

det(i + S) = l + tr& + ]Tdet (2) &( 2 )+ £ det (3 V 3 ) 

i<j i<j<k 

+ +det(S) , (All) 
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where e.g. det^ 3 -^ 3 ) denotes a 3 x 3 subdeterminant of b 
with the indices given by the summation variables. The 
comparison with the calculation of (e aA ) yields 



(e aA ) = det [i + (e" a - l)n ] 



(A12) 



This result also follows from Klich's trace formula 8 



Tr(e A e s ) 



dct(l 



(A13) 



(B4) 



with 



M 

7(2) = .9 6 °WE T 'l ' a )< ' a l > 

a=l 

where h r is the one particle Hamiltonian on the ring and 
g®(z) is the diagonal element of the resolvent of the semi- 
infinite chain at the boundary. For z = e ± iO it is given 
by 



where A and B are arbitrary one particle operators in 
Fock space and a and b are the corresponding operators 
in the Hilbert space of a single particle. The proof of this 
formula is nontrivial when the operators A and B do not 
commute. As Klich's elegant proof 8 - involves a result not 
widely known and it is not availble in a journal publi- 
cation, we presented the alternative proof of Eq. (|A12[) 
using Wick's theorem to make this paper selfcontained. 
In section II we use a slight generalization of Eq. (|A12|1 . 

Let C be a one particle operator in Fock space which 
commutes with Hq. With H Q = Hq — "fC we obtain 



(e aA e lC ) 



Tre aA e~ S ° Tre _J? o 



Tre 



-H' 



Tre 



det 



det(i 



det(l 



det [i + (e° 



l)no] 



(A14) 



APPENDIX B: PROOF OF EQ. |(77|) 

The scattering states can expressed via the full resol- 
vent g(z) = (z — h)^ 1 and the unperturbed resolvent 
go(z) = [z — ho) which are related by g — go + gvgo- 
From the definition in Eq. (|30[) one obtains for a! ^ a 

(m,a\k,a+) = (m, a\go(ek + iO)\l, a)r a (Bl) 
x (0,a\g(e k + iO)\0,a')T a ,(0,a'\k,a') . 

The resolvent matrix element of the semi-infinite chain is 
given by^ 3 - (m, a\g (e k + iO)|l,a) = - e ~ lkm . With Eqs. 
(f2~Tj) . (f2"2"| and ([57)1 the transmission probabilities follow 



h °(e±zO) = (eT*V / 4-e 2 )/2. 



(B5) 



In order to prove the relation between the transmission 
probabilities in Eq. (]TT[) we us a simple operator identity. 
Let a and b operators such that the inverses of a± b exist. 
The one has 



1 



1 



a — b a + b a — b 



1 / LN 1 1 / o 1 

( a + b )^TT-—^( a - b ) 



a + b a — b 



2^b^- 



a + b 
(B6) 



If one puts the unit operators (a ± b)(a ± on the 
opposite sides one obtains by comparison the operator 
relation 



1 b^ = ^b 1 



a — b a + b a + b a — b 



(B7) 



which we now use for operators acting in the subspace 
of the ring states. With a = el — h r — Re7(e + zO) and 
b = Im7(e + iO) = 71(e) the relation reads using Re7(e + 
iO) = Re7(e — iO) and Im7(e + iO) = — Im7(e — iO) 

g(e + iO) 7 i(e). 9 (e - iO) = g(e - iO) 7l (e) 5 (e + iO). (B8) 

We next take the expectation value in the state |0, M) 
and use the explicit form of j(z) in Eq. (|B4j) . The M = a 
term of the sum can be omitted on both sides of the 
equation as they are identical. Deviding by g°(e + iO) 
yields 



I* 



a < — a \2 



4r a 2 r Q 2 ,sin 2 fc|(0,a'| 5 (e fc + i0)|0,a)| 2 , (B2) 



Using the projection P r onto the states on the ring the 
full resolvent matrix elements in Eq. (|B2|) can be written 



(0,o'|<7(z)|0,o) = {0,a'\[zP r 



7 (z)]- 1 |0,a) (B3) 



AI-l 



M-l 



J2 r 2 |(0,A/| 5 (e + l 0)|0,a)| 2 = ^ r 2 |(0, a\g(e + i0)\0, M)f 



a=l 



Together with the explicit result for the transmission 
probalities in Eq. (|B2[) this completes the proof of the re- 
lation between the transmission probabilities in Eq. (|77|) . 
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